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ABSTRACT 



Sagittarius A*, the super- massive black hole at the center of the Milky Way, is surrounded by a small cluster of high 
velocity stars, known as the S-stars. We aim to constrain the amount and nature of stellar and dark mass associated 
with the cluster in the immediate vicinity of Sagittarius A*. We use near-infrared imaging to determine the _K's-band 
luminosity function of the S-star cluster members, and the distribution of the diffuse background emission and the stellar 
number density counts around the central black hole. This allows us to determine the stellar light and mass contribution 
expected from the faint members of the cluster. We then use post-Newtonian N-body techniques to investigate the effect 
of stellar perturbations on the motion of S2, as a means of detecting the number and masses of the perturbers. We find 
that the stellar mass derived from the i^s-band luminosity extrapolation is much smaller than the amount of mass that 
might be present considering the uncertainties in the orbital motion of the star S2. Also the amount of light from the 
fainter S-cluster members is below the amount of residual light at the position of the S-star cluster after removing the 
bright cluster members. If the distribution of stars and stellar remnants is strongly enough peaked near Sagittarius A*, 
observed changes in the orbital elements of S2 can be used to constrain both their masses and numbers. Based on 
simulations of the cluster of high velocity stars we find that at a wavelength of 2.2 /^m close to the confusion level for 
8 m class telescopes blend stars will occur (preferentially near the position of Sagittarius A*) that last for typically 3 
years before they dissolve due to proper motions. 

Key words. Galaxy: center - infrared: general - infrared: diffuse background - stars: luminosity function, mass function 
- stars: kinematics and dynamics - methods: numerical 



1. Introduction 

Using 8-10 m class telescopes, equipped with adaptive op- 
tics (AO) systems, at near-infrared (NIR) wavelengths has 
allowed us to identify and study the closest stars in the 
vicinity of the super-massive black hole (SMBH) at the 
center of our Milky Way. These stars, referred to as the 
S-star cluster, are located within the innermost arcsecond, 
orbiting the SMBH, Sagittarius A*(Sgr A*), on highly ec- 
centric and inclined orbits. Up till now, the trajectories 
of about 20 stars have been pr ecisely determined us ing 
NIR imaging and spectroscopy (jGillessen et al.l l2009al lbl) . 
This orbital information is used to determine the mass of 
the SMBH and can in principle be used to detect rela- 
tivistic effects and /or the mass distributio n of the central 
stellar cluster fe ubilar & Eckart 2001; Zu cker et al.ll2006t 
iMouawad etloBoOSl: iGillessen et al.ll2009alb 

One of the brightest members of that cluster is the 
star S2. It has the shortest observed orbital period of 



* Based on observations collected at the European 
Organisation for Astronomical Research in the Southern 
Hemisphere, Chile (ProgId: 073.B-0085) 



^15.9 years, and was the star used to precisely deter- 
mine the enclosed dark mass, and infer the existence of 
a ~4 mil lion solar mass SMBH, in our own Ga lactic cen- 
ter (GC; iSchodelet ai][200l iGhez et all [ 200l The first 
spectro scopic stu d ies of S2, bv lGhez et al.l ()2003f ) and later 
Eiscnha uer et al.l (pMH), revealed its rotational velocity to 
be that of an 08-BO young dwarf, w ith a mass of 1 5 M(?, 
and an age of less than 10^ yrs. Later. [Martins et al.l (|2008f) 
confined the spectral type of S2 to be a BO-2.5 V main- 
sequence star with a zero-age main-sequence (ZAMS) mass 
of 19.5 Mq. The fact that S2, along with most of the S- 
stars, is classified as typical solar neighborhood B2-9 V 
stars, indicates t hat they are yo ung, with ages between 6- 
400 Myr (Eiscnh auer et al.ll2005D . The combination of their 
age and the proximity to Sgr A* presents a challenge to star 
formation theories. It is still unclear how the S-stars were 
formed. Being generated locally requires that their forma- 
tion must have occurred through non-standard processes, 
like formation in at least one gaseous disk (Lockmann ct al] 
|2009| ) or via a n eccentricity instab ility of stellar disks 
around SMBHs (jMadigan et al.ll2009D . Alternatively, if they 
formed outside the central star cluster, about 0.3 par- 
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Fig. 1. Left: 17.5" x 17.5" NACO _R's-band mosaic of the central cluster zoomed in to the inner 1-2" region around Sgr A*. The 
inner region is indicated by a dashed (white) circle. Right: a) Map of the diffuse background light within a circle of 0.69" radius 
centered on the position of Sgr A*, shown here as a cross at the center. The projected orbit of the star S2 is over-plotted as an 
ellipse, b) The same map smoothed by convolution with a Gaussian beam of FWHM — 6 pixels. The contours levels are at 95%, 
90%, 80%, 70%, 60%, 50%, 40%, 30%, 20% and 10% of the maximum flux value for each image. 



sec c ore radius (e.g. iBuchholz et al.l 120091 : ISchodel et al.l 
|2007() , there are several models that describe how they ma 
have been brought in (e.g . JJansen & Milosavlicvic 200. 



Kim et al. 2004; Lev: 



i 



n et all 2005: Fujii ct al. 2009, 2010: 



Merritt et al.l l2009t iGould fc Quillen, ,2003; .Perets et al 



20071 20091). For a detai l ed de scription of these processes 
see 'Perets fc GualandrisI (|2010[ ). 

Stellar dynamics predict the formation of a cusp of 



stars at the center of a relaxed stellar cluster around 
a SMBH. This is manifested by an increase in the 
three dimensional stellar density of old stars and rem- 
nants towards the center with power-law slop es of 
T5 to 1.75 (iBahcah fc Wolj [19761 iMurphv et al.l [T99lt 



iLightman k Shapirolll977t [Alexander fc HoDmanll2009D . 

The steep power-law slope of 1.75 is reached in the case 
of a spherically symmetric single mass stellar distribution 
in equilibrium. For a cluster with differing mass compo- 
sition, mass segregation sets in, where the more massive 
stars sink towards the center, while the less massive ones 
remain less concentrate d. This leads to the s hallow den- 
sity distribution of 1.5 (jBahcah fc Wolj FT977I) . Later nu- 
merical simulations and analytical models confirmed these 
results (iFreit ag et all 120061: [Preto fc Amaro-Seoand [20Tot 



iHopman fc A lcxande: 



[| l2006l:IP 
a1 l2006bD . 



These steep density distri- 



butions were expected for the central cluster considering 
its age, which is comparable to the estimates of the two- 
body relaxation- t ime of 1-20 Gyr for the central parse c 
(| Alexander! I2005t lMerrit^l2010t iKocsis fc Tremainell201lD . 
However, observations of the projected stellar number den- 
sity, which can be related to the three dimensional density 
distribution, revealed that the cluster's radial profile can be 
fitted by two power-law slopes. The slope for the whole clus- 
ter outside a radius of ~ 6" (corresponding to 0.22 parsec) 
was found to be as steep as 1.8 ±0.1, while inside the break 
radius the slope was shal lower than expected and reached 
an ex ponent of 1.3 ± 0.1 (jGenzel et al.ll2003l : ISchodel et al.l 
l2007t) . These findings motivated the need to derive the den- 
sity profiles of the distinct stellar populations, given that re- 
cent star formation (6 Myr, Paumard ct al. 2006) at the GC 
gave birth to a large number of high-mass young stars that 
would be too young to reach an equilibrium state. Using 
adaptive optics and in termediate-band spectrophotometry 
iBuchholz et al.l (|2009f ) found the distribution of late-type 
stars (K giants and later) to be very flat and even showing 
a decline towards the Genter (for a radius of less that 6"), 
while the early-type stars (B2 main-sequence and earlier) 



fol low a steeper pr ofile. Similar results were obtained later 
bv lDo et al.l (|2009 ^ and'B artko et al.l (|2010[ ). 

These surprising findings required new models to ex- 
plain the depletion in the number of late-type giants in 
the central few arcseconds around the SMBH. Such at- 
tempts involved Smooth Particle Hydrodynamics (SPH) 
and Monte Carlo simulations which tried to account for 
the under density of giants by means of collisions with othe r 
stars and stellar remnants (|Dale et al.ll2009tlFreitaell2008[ ). 
Another explanation could be the disturbance of the cusp 
of stars after experiencing a minor merger event or an in- 
spiraling of an intermediate-mass black hole, which then 
would lead to deviations from equilibriur n; hence causing a 
shallower power-law profile of the cusp (jBaumgardt et al.l 
|2006| ) . iMerrittI (|2010D explains the observations by the evo- 
lut ion of a parsec-s c ale in itial core model. 

iMouawad et al.l (|2005[) presented the first efforts to de- 
termine the amount of extended mass in the vicinity of 
the SMBH allowing for non-Keplerian orbits. Using posi- 
tional and radial velocity data of the star S2, and leav- 
ing the position of Sgr A* as a free input parameter, they 
provide, for the first time, a rigid upper limit on the pres- 
ence of a possible extended dark mass component around 
Sgr A*. Gonsidering only the fraction of the cusp mass 
Msi 2,^„ that ma. y be within the apo-center of the S2 orbit, 
iMouawad et al.l 12005 ) find Ms2,p„/(Msmbh + Mga^^J < 
0.05 as an upper limit. This number is consistent with 
more recent investigations of the problem (jGillessen et al.l 
l2009b[ ). Due to mass segregation, a large extended mass 
in the immediate vicinity of Sgr A*, if present, is unlikely 
to be dominated in mass of sub-solar mass constituents. It 
could well be explained by a cluster of high M/L stellar 
remnants, which may form a stable configuration. 

From the observational point of view, several at- 
tempts have been made recently to tackle the missing 
cusp problem. iSazonov et al.l ()201lD proposed that the de- 
tected 1" sized thermal X -ray emission close to Sgr A* 
(|Baganoff et al.ll200ll I2003D can be explained by the tidal 
spin-ups of several thousand late-type main-sequence stars 
(MS). They use the Chandra X-ray data to infer an up- 
per limit on the density of these low-mass main-sequence 
stars. Furthermore, using Hub ble Space Telescope (HST) 
data, lYusef-Zadeh et all (J2012t ) derived a stellar mass pro- 
file, from the diffuse light profile in the region < 1" around 
Sgr A*, and by that they explained the diffuse light to be 
dominated by a cusp of faint KO dwarfs. 
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Up to now, the true distribution of the Nuclear Star 
Cluster, especially the S-stars, is yet to be determined. No 
investigations have confirmed or ruled out the existence of 
a cusp of relaxed stars and stellar remnants around Sgr A*, 
as predicted by theory. An excellent dataset to investigate 
the stellar content of the central arcsecond around Sgr A* is 
the NIR i^s-band (2.2 fj,m) data (see Figure [TJ we used in 
(jSabha et al.l 120101 hereafter NSIO). In that case we sub- 
tracted the stellar light contribution to the flux density 
measured at the position of Sgr A*. The aim of this work 
is then to analyze the resulting image of the diffuse NIR 
background emission close to the SMBH. This emission 
is bel ieved to trace the accumulative light of un resolved 
stars (jSchodel et al.ll2007trYu"sef-Zadeh et al.ll2012[ ). We ex- 
plain the background light by extrapolating the iiTs-band 
luminosity function (KLF) of the innermost (1-2", corre- 
sponding to 0.05 parsecs for a distance of 8 kpc to the GC) 
members of the S-star cluster to fainter iCs-magnitudes. We 
compare the cumulative light and mass of these fainter stars 
to the limits imposed by observations. We then extend our 
analysis to explore the possible nature of this background 
light by testing its effect on the observed orbit of the star 
S2. Furthermore, we simulate the distribution of the unre- 
solved faint stars {Kg > 18) and their combined light to 
produce line-of-sight clusterings that have a compact, close 
to stellar, appearance. 

The paper is structured as follows: Section [2] deals with 
a brief description of the observation and data reduction. 
We describe in Section|3]the method used (? 13.1] -? [53)1 and 
discuss the different observational limits (§ 13. 4p employed 
to test our analysis. Exploring the possible contributors to 
the dark mass within the orbit of S2 is done in Section^ In 
Section[S]we give the results obtained by simulating the dis- 
tribution of faint stars and the possibility of producing line 
of sight clusterings that look like compact stellar objects. 
We summarize and discuss the implications of our results 
in Section[6l We adopt throughout this paper S(i?) ex R^^ 
as the definition for the projected density distribution of 
the background light , with R being the projected radius 
and r the corresponding power-law index. 

2. Observations and data reduction 

The observations and data reduction have been described in 
NSIO. In summary: The near-infrared (NIR) observations 
have been conducted at the Very Large Telescope (VLT) 
of the European Southern Observatory (ESO) on Paranal, 
Chile. The data were obtained with YEPUN, using the 
adaptive optics (AO) module NAOS and the NIR cam- 
era/spectrometer CONIC A (briefly "NACO"). The data 
were taken in the Xg-band (2.2 /itm) on the night of 23 
September 2004, and is one of the best available where 
Sgr A* is in a quiet state. The flux densities were mea- 
sured by aperture photometry with circular apertures of 
66 mas radius. They were corrected for extinction, us- 
ing Afc„ = 2 .46 d erived for the inner arcsecond from 
ISchodel et al.l ()2010[ ). Possible uncertainties in the extinc- 
tion of a few tenths of a magnitude do not influence the gen- 
eral results obtained in this paper. The flux density calibra- 
tion was carried out using zero points for the correspond- 
ing camera setup and a comparison to kno wn j^g-band flux 
dens ities of IRS16C, I RS16NE (from Schodel et al.l 120101: 
also Blum et al.l 19961 ) and to a number of the S-stars 
(|Witzel et al.ll2012D . 



3. The central few tenths of parsecs 

In NSIO we gave a stringent upper limit on the emission 
from the central black hole in the presence of the sur- 
rounding S-star cluster. For that purpose, three indepen- 
dent methods were used to remove or strongly suppress 
the flux density contributions of these stars, in the central 
~ 2", in order to measure the flux density at the posi- 
tion of Sgr A*. All three methods provided comparable re- 
sults, and allowed a clear determination of the stellar light 
background at the center of the Milky Way, against which 
Sgr A* has to be detected. The three methods, linear ex- 
traction of the extended flux density, automatic and itera- 
tive point spread function (PSF) subtraction were carried 
out assuming that the extracted PSF in the central few 
arcseconds of the image is uniform . Investigations of larger 
images (e.g. iBuchholz et al.ll2009t ) show that on scales of a 
few arcseconds the constant PSF assumption is valid, while 
for fields > 10" the PSF variations have to be taken into 
account. 

Figured] is a map of the 51 stars adopted from the list 
in Table 3 of NSIO. The stars are plotted relative to the 
position of Sgr A*. The surface number density of these 
detected stars, within a radial distance of about 0.5" from 
Sgr A*, is 68±8 arcsec"^, with the uncertainty correspond- 
ing to the square-root of that value. This value agrees with 
the central num ber density of 60 ± 10 arcsec"^ given by 
iDo et al.l (|2009D . Extrapolating the KLF allows us to test if 
the observed diffuse light across the central S-star cluster, 
or the amount of unaccounted dark mass, can be explained 
by stars. 
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Fig. 2. Map of the 51 stars listed in Table 3 from NSIO. The 
color of each star indicates its A'a-magnitude. The size of each 
symbol is proportional to the flux of the corresponding star. The 
position of Sgr A* is indicated as a cross at the center. 
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Fig. 3. KLF histogram of the stars detected in the central field, 
derived from the 23 September 2004 data. The dashed line in- 
dicates the linear fit of the KLF slope of a = 0.18 ± 0.07. The 
vertical dotted line (red) represents the current detection limit 
for faint ifs-magnitudes. 



3.1. KLF of the S-star cluster 

Figure [3] shows the KLF histogram derived for the stars 
detected in the central field, (Figure [5]). We improve the 
KLF derivation by choosing a fixed number of bins that 
allows for about 10 sources per bin while providing a suf- 
ficient number of points to allow for a clear linear fit. The 
Red Clump (RC) /Horizontal Branch (HB) stars, around 
Kg « 14. 5, are in one bin, so the RC/HB bump is visi- 
ble there (jSchodel et al.l 120071 ). For estimating the uncer- 
tainty, we randomized the start of the first bin in an in- 
terval between Kg = 13.0 to 14.2 and repeated the his- 
togram calculation 10^ times. The number of sources in 
each bin was then determined by taking the average of all 
iterations and the uncertainties were subsequently derived 
from the standard deviation. We derive a least-square linear 
slope oi dlog{N)/d{Ks) — a = 0.18±0.07, which compares 
well with the KLF slope of 0.3 ± 0.1 derived in NSIO and 
also with the KLF slope of 0.21 ± .02 fo und for the inner 
field {R < 6") bv iBuchholz et al.l (|2009h . For the magni- 
tudes up to Kg = 17.50 within the central 0.69" radius, 
we detect no significant deviation from a straight power- 
law. This implies that the completeness is high and can be 
compared to t he ^7 0% value derived for mag^ = 17 by 
ISchodel et al.l ([2003) where the authors introduced artifi- 
cial stars into their NIR image and attempted re-detecting 
them. However, for K^ — 17.50 to 18.25 the stellar counts 
drop quickly to about 20% of the value expected from the 
straight power-law line; hence the last i^Tg-bin is excluded 

fro m the linear fit. 

iMaiz Apellaniz fc Ubedal ([200^ propose an alternative 
way of binning when dealing with stellar luminosity and ini- 
tial mass functions (IMF). Their method is based on choos- 
ing variable sized bins with a constant number of stars in 
each bin. They find that variable sized binning introduces 
bias- free estimations that are independent from the number 



of stars per bin. Their method is applicable to small samples 
of stars. We apply their method to our KLF calculation and 
get d{log{N)/6Ks)/d{Ks) = 0.12±0.09, consistent with our 
fixed sized binning method. 

3.2. The diffuse NIR background 

The methods we used in NSIO to correct for the flux density 
contribution of the stars in the central 2" have revealed a 
faint extended emission around Sgr A* (NSIO Figures 3b, 
4b and 5). We detected ~ 1.3 mJy (obtained by correct- 
ing the ~ 2 mJy we quote in NSIO for the Ak^ = 2.46 
we use here) at the center of the S-star cluster. With a ra- 
dius of 1" (about twice the FWHM of the S-star cluster) 
for the Point Spread Function (PSF) used for the subtrac- 
tion, we showed that a misplacement of the PSF for about 
only five stars, located within one FWHM of Sgr A*, would 
contribute significantly to the measured flux at the center. 
For a median brightness of about 1.3 mJy for these stars, a 
1 pixel ~ 13 mas positional shift of each of these stars to- 
wards Sgr A* would be required to explain all the detected 
^ 1.3 mJy at the center i.e. 0.26 mJy from each star. In 
ISabha et al.l (|201lh we showed that a displacement larger 
than a few tenths of a pixel would result in a clear and 
identifiable characteristic plus/minus pattern in the resid- 
ual flux distribution along the shift direction. For a maxi- 
mum positional uncertainty of 1 pixel, we showed that the 
independent shifts of the five stars can be approximated 
by a single star experiencing five shifts in a random walk 
pattern. This resulted in calculating a total maximum con- 
tribution of 0.26 mJy from all the five stars to the center, 
which translates to about 20-30% of the flux density. Thus, 
more than two thirds of the extended emission detected to- 
wards Sgr A* could be due to faint stars, at or beyond the 
completeness limit reached in the KLF, and associated with 
the ~ 0.5-1" diameter S-star cluster. 

The diffuse background emission we detected (see 
Figure [IJl) could be compared to the projected distribution 
of stars S(i?) oc R~^ , with R being the projected radius. 
We found that the distribution of the azimuthally averaged 
residual diffuse background emission, centered on the po- 
sition of Sgr A*, not to be uniform but in fact decreases 
gently as a function of radius (see Figure 7 in NSIO) with 
a power-law index Fdiftusc = 0.20 ± 0.05. In this investiga- 
tion we re-calculate the azimuthally averaged background 
light from the iterative PSF subtracted image alone. The 
azimuthally averaged background light is plotted as a func- 
tion of projected radius from Sgr A* in FigureU) In this new 
calculation we find the power-law index to have a value of 
FdiffusG = 0.14 ± 0.07. Both results are consistent with re- 
cent investigations concerning the distribution of number 
density counts of the stellar populations in the central arc- 
seconds, derived from im aging VLT and Keck data. For 
the ce ntral few arcsecond s lBuchholz et al.l (|2009[ ). IDo et al.l 
(I2OO9I) and iBartko et all ([20101) find a F ~ 1.5 ± 0.2 for 
the young stars, but an even shallower distribution for the 
late-type (old) stars with F ~ 0.2 ± 0.1. A detailed discus- 
sion concerning the different popula tions and their distribu- 
tion is given inlGenzel et al.l (l2003l): ISchodel et al.l (I2007D: 



IBuchholz eral! (|2009[ ): iDo et all (|2009t ) and iBartko et all 
(I2010D . 

The small value we obtain for the projected diffuse 
light exponent Fdiffuse and the high degree of completeness 
reached around Kg = 17.5, makes this data set well suited 
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for analyzing the diffuse background light. Especially in in- 
vestigating the role of much fainter stars, beyond the com- 
pleteness limit, in the observed power-law behavior of the 
background. 
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Fig. 4. Azimuthal average of the diffuse background emission as 
derived from manual PSF subtracted 23 September 2004 image. 
The squares (mean flux and Icr uncertainty per pixel) have been 
calculated in annuli of 39.8 mas (3 pixels) width. The black 
dashed line marks a fit to the data points with an exponential 
decrease of 0.14. 



3.3. Extrapolating the KLF of the S-star cluster 

Motivated by the power-law behavior of the diffuse back- 
ground emission and assuming that the drop in the KLF 
counts at magnitude ~ 18 is caused only by the fact that we 
have reached the detection limit, we extrapolate the KLF to 
fainter magnitudes in order to investigate how these faint 
stars contribute to the background light. The true shape 
of the luminosity function for i^'g-magnitudes below the 
completeness limit of ~ 17.5 has yet to be determined. 
Investigations into the IMF of the S-cluster have shown 
that it can be fitted with a standard Salpeter/Kroupa IMF 
of dN /dm ex m^^'^ and continuous sta r formation h i storie s 
with moderate ages (below 60 Myr, iBartko et al.l |2010|) . 
Here, we estimate an upper limit on the stellar light by as- 
suming that the KLF exhibits the same behavior observed 
for brighter magnitudes without suffering a break in the 
slope toward the fainter end. 

We use the KLF slope we obtained for the innermost 
central region, 0.18 ± 0.07 (Figure [3]) and extrapolate it 
over five magnitudes bins to Ks ~ 25. The Xg-magnitude 
bins between 18-25 (translating to stellar masses in the 
range of ~ 1.68 to 0.34 Mq) correspond to the bright- 
ness of the expected main-sequence stars (luminosity class 
V) which are likely to be present in the central cluster. 
However, we assume tha t due to mass se gregation effects 
in th e Galactic nucleus (jBahcall fc Wolj 1976: Alexander: 
|2005[ ). driven by dynamical friction ( Chandrasekharn 19431) 



between stars, the heavier objects sink towards the center 
while the lighter objects move out. Their volume density 
will be significantly re duced and they may e ven be expelled 
from the very center. iFreitag et al.l (J2006l) show that the 



main-sequence stars begin to be expelled outward by the 
cusp of stcUar-mass black holes (SBH) after a few Gyrs, 
just shorter than the presumed age of the stellar cluster 
at about 10 Gyrs. While the reservoir of lower mass stars 
may be replenished by the most recent - possibly still on- 
going - star formatio n episode about 6 million years ago 
(jPaumard et al.ll2006[ ). we assume that stars well below our 
low mass limit of ^ 0.34 Mq with Kg-hand brightnesses 
around K^ — 25 are affected by depletion. 

Figure [5] shows the KLF slope of a = 0.18 and the 
upper limit imposed by the uncertainty in the fit (a — 0.25) 
plotted as dashed and dash-dotted lines, respectively. The 
extrapolated Xg-bins are shown as hollow circles. We adopt 
a Monte Carlo approach for calculating the number of stars 
A'' from the KLF, taking into account the uncertainty in the 
slope. After 10^ trials we find as a result for each bin, the 
median number N and median deviation dN . 
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Fig. 5. Extrapolation of the KLF power-law fit. The KLF slope 
of Q = 0.18 and the upper limit imposed by the uncertainty in 
the fit (a = 0.25) are plotted as dashed and dash-dotted lines, 
respectively. The black filled circles represent the data while the 
hollow circles represent new points based on the extrapolated 
KLF slope. The approximate location of the detection limit is 
indicated by the vertical dotted/red line. 



Using the extrapolated -Ks-magnitudes, the correspond- 
ing flux densities are calculated using the following relation 



/n 



= /S2 X 10-"-4 (^" 



-Ks2) 



(1) 



where /now star and K^^^w star are the flux density and 
iiTs-magnitude for each new star in the extrapolation. The 
flux and magnitude for the star S2 were adopted from NSIO, 
Table 3, and corrected for the extinction value we use here 
(see § [2]). The new values are /s2 = 14.73 mJy and Kv,2 = 
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14.1. The accumulative flux density for each Ks-hin /bii 
obtained via 



IS 



/bin — In 



XNr, 



(2) 



The number of stars per bin A^ncw star is randomly picked 
from the interval between [iVncw star — dN^cw star] and 
[.^new star + dN^cw star]- In 10^ trials the accumulative flux 
per bin and its uncertainty are determined as the median 
and median deviation of the randomly drawn fluxes /bin- 
We then add up all the accumulative flux densities for all 
the new Ks-hand bins and obtain the integrated brightness 
of the extrapolated part of the S-star cluster, 



25 



i^Extra stars = ^ /bin = (25.72 ± 14.31) mjy. 



(3) 



K,~18 



We assume that the faint, undetectable stars follow the 
distribution of the azimuthally averaged background light, 
as shown in Figure ID Thus, the light from the faint stars 
that we introduced in the 0.69" radius region can be com- 
pared to the measured background light from our data 
for the same region. This is achieved by using the total 
flux density Fsxtra stars to derive the peak light density 
(^Extra Stars) that would bc mcasurcd inside one resolution 
element of 0.033" radius centered on the position of Sgr A*, 
using the following relation: 



Extra Stars 



/(r, (f))rdrd(p 



2ttL 



Extra Stars 



0.690 



0.033 



„i~r 



dr, 



with F = Fdiffuso = 0.14 (see § E^]). The peak hght den- 
sity for the extra stars is then /Extra stars — (15.24 ± 
8.48) mJy arcsec"^. To compare the hght caused by the 
extra stars with the measured background emission, we 
plot the stellar light density caused by our new stars with 
the azimuthally averaged measured light density of the 
background (Figure [5]) . For illustration purposes we nor- 
malize the observed peak stellar light to the measured 
background value within the central resolution element, 
^Baciiground = (254.30 ± 58.45) mJy arcsec"^. It is clear 
that the peak light introduced by the new faint stars, as 
calculated from the extrapolation of the 0.18 ± 0.07 KLF 
slope, is very small and below that of the background. The 
dotted line (black circles) represents the background light 
while the dashed line (blue squares) corresponds to the ex- 
tra stellar light. The upper limit of the extrapolated extra 
stellar light contribution is presented as a dashed line with 
no symbols. The figure shows that the upper limit of the 
extrapolated light contribution of the S-star cluster is lower 
than 15% of the measured background light. 



3.4. Observational limits on the stellar light and mass 

Our analysis shows that if there was a population of very 
faint stars, following the extrapolated Kg-hand luminosity 
function and central cluster profile obtained for the brighter 
stellar population (less than Kg — 18), the additional stellar 
light and mass lie well below the limits given by observa- 
tional data. See following sections and Figures [7] and \E\ 
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Fig. 6. Relative azimuthally averaged hght density, for the 
background light taken from the observations and the extra stel- 
lar light calculated from the extrapolation, plotted as a function 
of distance from Sgr A*. They are represented by a dotted line 
with circles (black) and dashed hnes with squares (blue), respec- 
tively. The stellar light density is normalized to the peak light 
density of the background at the central resolution element. The 
upper limit of the extrapolated extra stellar light is shown as the 
blue dashed line with no symbols. 



(4) 3.4.1. Limits on the stellar light 



Following the previous calculations and the result displayed 
in Figure El we perform our analysis for a range of KLF 
slopes in order to test if the observed background light can 
be solely obtained by the emission of faint stars. The range 
of KLF slopes we use is based on the values and uncer- 
tainty estimates of the following published KLF slopes for 
the central 2": 0.13 ±0.02 (Buchholz et al. 2009, early-type 
stars), 0.27 ± 0.03 (iBuchholz et al.ll2009l . late-tvpe stars), 
0.21 ± 0.02 (jBuchholz et al.ll2009l all stars) and 0.30 ±0.1 
(NSIO), in addition to the improved newly fitted slope of 
the KLF in this work 0.18 ± 0.07. 

We extrapolate each KLF slope to a /sTg-magnitude of 
~ 25. The peak light density (/Extra stars) is calculated us- 
ing Equation @ . The peak light density of the extra stars is 
plotted for the extrapolated KLF slopes in the range of 0.11 
to 0.40 in Figure [7] . The limit imposed by the peak light 
density of the measured background light (Figure [J) is plot- 
ted as a horizontal dashed line (blue). In ad dition, the KLF 
slopes derived in this work and by NSIO andl Buchholz et aTl 
(200il) are plotted as purple, yellow and green data points, 
respectively. Figure [7] clearly shows that almost all of the 
KLF slopes result in a peak light density below the ob- 
served limit, except for very high slopes > 0.37 which are 
not in agreement with the observations. 

3.4.2. Limits on the stellar mass 

Using the same range of KLF slopes, we estimate the mass 
that would be introduced to the central region as a result 
of the KLF extrapolation. We obtain the stellar mass cor- 
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Fig. 7. Estimated pealc light dens ity from stars d e rived from 
for different KLF slopes. Slopes of Buchholz et al.l (|2009l ') (for 
different stellar populations) are shown in green, NSIO in yellow 
and the KLF slope derived here, in § 13.11 in purple. A limit 
imposed by the measured peak light density from the measured 
background light is plotted as a horizontal dashed line (blue). 



Fig. 8. Estimated stellar ma ss from the add e d sta rs for dif- 
ferent KLF slopes. Slopes of iBuchholz et all (120091 ) (for dif- 
ferent stellar populations) are shown in green, NSIO in yellow 
and our fitted slope in purple. A limit imposed by the enclosed 
mass within the S2 orbit is plotte d as horizontal dotted (gray) 
and dash-dotted (red ) lines from iGillessen et al.l 12009b) and 
iMouawad et al.l (|2005r ). respectively. 



responding to the extrapolated i^Tg-bins by calculating their 
luminosity via 



Lk^ - lO-0-MM^.-MeK^)j^^ 



QK,, 



(5) 



where, Lk, and Mk^ are the luminosity of a star and its ab- 
solute magnitude in i^Tg-band, respectively. Lqk, & Mqk^ 
are the Kg luminosity and absolute magnitude of the Sun. 
Then, the mass for each Kg magnitude is calculated using 



m = (LkJ(V4) 



(6) 



from iDurid (|2004 ): ISalaris fc Cassisil ^M)- For example, 
a ifg-magnitude around 20 corresponds to 1 Mq main- 
sequence stars of FOV, GOV, K5V spectral types. 

In Figure \8\ we show the estimated extra mass for all 
the KLF slopes in units of solar mass. The figure also 
shows, dash-dotted/red line , the upper limit for an ex- 
tended ma ss enclosed by the orb it of the star S2, cal- 
culated by IMouawad et al.l (|2005D . where they use non- 
Keplerian fitting of the orbit to derive the upper limits, 
assuming that the composition of the dark mass is sources 
with M/L ^ 2. The dotted/gr ay line represeii t s the ti ghter 
upper limit obtained later by IGillessen et al.l (J2009bD who 
derive the mass using recent orbital parameters of S2. They 
assum e that the extende d mass consists of stellar black 
holes (iFreitag et alJ 20061) wi t h a in ass o f 10 Mp) using est i- 
mations from iTimmes et al.l (|l996l ) and lAlexanden (J2007f) . 
It can be concluded from the figure that the introduced 
stellar mass, within a radius of ~ 0.69", lies well below the 
upper limits imposed by the S2 orb it with a semi-major axis 
of- 0.123" (jGillessen et al.ll2009bD . See Figured] (right) for 
a comparison of the sizes of the two regions. 



4. Dynamical probes of the distributed mass 

If the gravitational force near Sgr A* includes contributions 
from bodies other than the SMBH, the orbits of test stars, 
including S2, will deviate from Keplerian ellipses. These de- 
viations can be used to constrain the amoiint of distributed 
mass near Sgr A* (jMouawad et al.l 120051 : IGillessen et al.l 
l2009b[ ). But they can also be used to constrain the "gran- 
ularity" of the perturbing potential, since the nature and 
magnitude of the orbital deviations depend both on the to- 
tal mass of the perturbing stars, and on their individual 
masses. 

Investigations of a single scattering event were explored 
by Gualandris et al. (2010) using high-accuracy N-body 
simulations and orbital fitting techniques. They found that 
an IMBH more massive than 10'^ Mq, with a distance com- 
parable to that of the S-stars, will cause perturbations of 
the orbit of S2 that can be observed after the next peri- 
bothrorQ passage of S2. Here we examine the effect many 
scatterers (i.e. smaller masses for the scatterers but shorter 
impact parameters) will have on the trajectory of the star 
S2 as it orbits. Around Sgr A*, the stars and scatterers are 
moving in a potential well that is dominated by the mass of 
the central SMBH. In this case the encounters are of a cor- 
related nature and hence cannot be considered as random 
events. 

An important deviation from Keplerian motion occurs 
as a result of relativistic corrections to the equations of 



^ Peri- or apobothron is the term used for peri- or apoapsis for 
an elliptical orbit with a black hole present at the appropriate 
focus. 
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motion, which to lowest order predict an advance of the 
argument of peribothron, w, each orbital period of 



... _ 6ttGM, 



= 5.0 mpc and e = 0.88 for the semi-major axis 



Setting a = 5.0 mpc and e = 0.88 for the semi-maji 

i eccentricity of S2, respectively, and assuming M, — 
X 10^ Mq, ' 

(8) 



and 

4.0 X 10^ Mq 

(Aw)gr « 10.8'. 

The relativistic precession is prograde, and leaves the ori- 
entation of the orbital plane unchanged. 

The argument of peribothron also experiences an ad- 
vance each period due to the spherically-symmetric compo- 
nent of the distributed mass. The amplitude of this "mass 
precession" is 



(Aw)M = "27rGM(e,7)Vr 



M^r < a) 



M, 



(9) 



Here, M^ is the distributed mass within a radius r = a^ and 
Gm is a dimensionless factor of order unity that depends 
on e and on t he power-law index of the density, p ex r~'^ 
(|Merrittll201^ . In the special case 7 = 2, 



G 



M 



(l + Vl-e^) 



0.68 for S2 



so that 



(Ao.) 



M 



-1.0' 



M^{r < a) 
103 M, 







(10) 



(11) 



Mass precession is retrograde, i.e., opposite in sense to the 
relativistic precession. 

Since the contribution of relativity to the periboth- 
ron advance is determined uniquely by a and e, which 
are known, a measured Aw can be used to constrain the 
mass enclosed within S2's orbit, by subtracting (Au;)gr 
and comparing the result with Equation (|lip . So far, this 
technique has yielded o nly upper limits on M^, of ^ 10^^ M, 
(|Gillessenet al.ll2009bD . 

The granularity of the distributed mass makes itself 
felt via the phenomenon of "resonant relaxation" (RR ) 
(|Rauch fc Tremaind [Tool iHopman fc Alexandeij [2006aD . 
On the time scales of interest here, orbits near Sgr A* re- 
main nearly fixed in their orientations, and the perturbing 
effect of each field star on the motion of a test star (e.g. S2) 
can be approximated as a torque that is fixed in time, and 
proportional to m, the mass of the field star. The net effect 
of the torques from N field stars is to change the angular 
momentum, L, of S2's orbit according to 



ic M. P 



(12) 



where Lc — y/GM^a is the angular momentum of a circu- 
lar orbit having the same semi-major axis as that of the 
test star. (Equation [T^] describes "coherent resonant relax- 
ation" ; on time scales much longer than orbital periods, 
"incoherent" resonant relaxation causes changes that in- 
crease as ~ -y/At.) The normalizing factor K is difficult to 
compute from firs t principles but should be of order unity 
(jEilon et al.l l2009l) . Changes in L imply changes in both 
the eccentricity, e, of S2's orbit, as well as changes in its 



orbital plane. The latter can be described in a coordinate- 
independent way via the angle A9, where 



(7) cos(A6') 



Li • L2 
L1L2 



(13) 



and {Li, L2} are the values of L at two times separated by 
At. If we set At equal to the orbital period of the test star, 
the changes in its orbital elements due to RR are expected 
to be 



|Ae|RR « K,Vn^ 



M,' 



(A0)rr « 27TKtVN 



M.' 



(14) 
(15) 



where N is the number of stars having a-values similar to, 
or less than, that of the test star and {K^, Kt} are constants 
which may depend on the properties of the field-star orbits. 

Because the changes in S2's orbit due to RR scale differ- 
ently with m and N than the changes due to the smoothly- 
distributed mass, both the number and mass of the per- 
turbing objects within S2's orbit can in principle be inde- 
pendently constrained. For instance, one could determine 
Af* — mN from Equations d?]) and ([TT|) and a measured 
Aoj, then compute myN by measuring changes in e or 
and comparing with Equations ([H]) or p3|) . 

We tested the feasibility of this idea using numerical in- 
tegrations. The models and meth ods were similar to those 
described in iMerritt et al.l (|2010[ ) . The N field stars were 
selected from a density profile n{r) ex r~^, with semi-major 
axes extending to Omax — 8 mpc. Initial conditions assumed 
isotropy in the velocity distribution. Two values for the field 
star masses were considered: m = 10 Af© and m = 50 Mq. 
One of the A^-body particles was assigned the observed 
mass and orbital elements of S2; this particle was begun at 
apobothron, and the integrations extended for one complete 
period of S2's orbit. Each of the A^ field-star orbits were in- 
tegrated as well, and the integrator included the mutual 
forces between stars, as well as post-Newtonian corrections 
to the equations of motion. The quantities Aw, Ae etc. for 
the S2 particle were computed by applying standard formu- 
lae to (r, v) at the start and end of each integration. 100 
random realizations of each initial model were integrated, 
allowing both the mean values of the changes, and their 
variance, to be computed. 

Figures [9l and ITOk show changes in w for S2. The me- 
dian change is well predicted by Equation pT|) . However 
there is a substantial variance. We identify at least two 
sources for this variance. (1) The number of stars inside 
S2's orbit difi'ers from model to model by ~ vN , resulting 
in corresponding changes to the enclosed mass, and hence 
to the precession rate as given by Equation (|TT|) . (2) When 
N is finite, the same torques that drive resonant relaxation 
also imply a change in the field star's rate of peribothron 
advance as compared with Equation (jlip . which assumes 
no tangential forces. While the dispersion scales roughly as 
•\/iV, as evident in Figure [HI the fractional change in Auj 
due to this effect scales as -- 1/a//V (jMerritt et al.ll2010[ ). 
Additional variance might arise from close encounters be- 
tween field stars and S2, and from the fact that the mass 
within S2's orbit is changing over the course of the integra- 
tion due to the orbital motion of each field star. 

Whereas the (average) value of Aw depends only on the 
mass within S2's orbit, the changes in e and 9 depend also 
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Fig. 9. Histograms of the predicted change in S2's argument 
of peribothron, uj, over the course one orbital period (~ 16 yr). 
The shift due to relativity, (Aa;)GR ~ 11', has been subtracted 
from the total; what remains is due to Newtonian perturbations 
from the field stars. Each histogram was constructed from inte- 
grations of 100 random realizations of the same initial model, 
with field-star mass m = 10 Mq , and four different values of the 
total number: iV = 200 (solid/black); N = 100 (dotted/red); 
A'' = 50 (dashed/blue); and A^ = 25 (dot-dashed/green). The 
average value of the peribothron shift increases with increasing 
Nm, as predicted by Equation pip . The reasons for the spread 
in Auj values are discussed in the text. 



on m^ as shown in Figures llOb and c. The lines in those 
figures are Equations ((T4)) and (fTS)) . with 



Ke = 1.4, Kt = 1.0. 



(16) 



(We have defined N in Equations ([TH) and ([15]) as the num- 
ber of field stars inside a radius of 9.4 mpc, the apobothron 
of S2.) For a given value of the enclosed mass, M^ = Nm, 
Figure [in] shows that the changes in e and 6 indeed scale as 
^ 1/VN or as ^ \/m, as predicted by Equations dTi]) and 

We can use these results to estimate the changes in 
oj, e and 9 expected for S2, based on theoretical models 
of the distribution of stars and stell ar remnants at the 
GC. In dynamically e volved models ('Freitag et al.l l2006t 
iHopman fc Alexander! l2006bt) , the total distributed mass 
within S2's apobothron, r « 10 mpc, is predicted to be ^ 
a few times 10'^ Mq. About half of this mass is in the form 
of main-sequence stars and half in stellar-mass black holes, 
with a total number N sa lO'^. When there are two mass 
groups, expressions like Equations p^ and (|15l) generalize 
to 



|Ae|RR = K. 
(A^)r,r, = 27:Kt 



mi^/Ni + 'm2^/N2 



M, 
mi\/Ni + m2\/N2 



M, 



(17) 
(18) 



assummg 

mi = 1 Mq, m2 = 10 Mq, A^i = 10^ N2 ^ 150 (19) 
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Fig. 10. Average values of the changes in to, e and 9 for S2 
over one orbital period ('^ 16 yr) in the A-body integrations. 
Filled circles are from integrations with tti = 50 Mq and open 
circles are for m = 10 Mq ; the number of field stars was 
A = {25, 50, 100, 200} for both values of m. The abscissa is the 
distributed mass within S2's apobothron, at r « 9.4 mpc. In 
each frame, the points are median values from the 100 A-body 
integrations, and the error bars extend from the 20th to the 80th 
percentile of the distribution, a) Changes in the argument of 
peribothron. The contribution from relativity. Equation (O, has 
been subtracted. The solid line is Equation (|ll|l . b) Changes in 
the eccentricity. Solid and dashed lines are Equation (|14|l . with 
m = 50 Mq and m — 10 Mq respectively and with Ke = 1.4. 
c) The angle between initial and final values of L for S2. Solid 
and dashed lines are Equations (|15|l with Kt = 1.0. 



(iHopman fc Alexandeij|2006bl ) we find 

|Ae|RR « 5.4 X 10-^ 
(A0)rr « 0'.8, 
(Aw)m « -2.5'. 



(20) 
(21) 
(22) 



For obtaining the dispersion in the value of Equation ([^^ , 
we scaled the dispersion given in Figure [TOk for the sin- 
gle population case, A^ = 50, m = 50 AIq of the same 
total extended mass, to the two populations case we are 
investigating here. The dispersion obtained from the sim- 
ulations is ~ 4'. We scale it using the relation Auj/\/N in 
order to account for the SBH and MS populations, indepen- 
dently. The dispersion for the new configuration then be- 
comes ~ 1.43', lower than the single population case. This 
is attributed to the fact that the number of main-sequence 
stars is much larger than the stellar-mass black holes, hence 
they lower the dispersion in the total Newtonian periboth- 
ron shift (Aa;)M- 
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Considering a higher value for the enclosed mass M^ — 
10^ Mq while keeping the same mass scales and abundance 
ratios of the scattering objects, 

mi = 1 Mq, ma = 10 Mq, 7Vi = 4000, iVs = 600 (23) 



one gets changes of 



|Ae|RR 
(A0)rr 

(Aw)m 



1.1 X 10" 

1.7', 

-10'. 



(24) 
(25) 
(26) 



|Ae|RR s 


a 3.4 X 10~^ 


(A0)rr s 


. 5.2', 


(Aw)m - 


s -100'. 



The dispersion in Equation (j26p can be compared, as 
we did before, to the case considered in the simulations 
(A^ = 200,771 = 50 M0) by scaling the ^ 8' dispersion 
(Figure ITUk) to become ~ 2.86' for the two mass popula- 
tion. 

Repeating the same analysis as before to the Af* = 
10^ Mq gives the following numbers for the stellar black 
holes and low-mass stars 

mi = 1 Mq, mz = 10 M0, iVi = 40000, N2 = 6000 (27) 

that result in 

(28) 
(29) 
(30) 

Similar to the above cases, the dispersion in Equation (j30p 
can be compared to the single mass case by scaling the 
^ 25.3' dispersion to become ~ 9.1' for the two mass 
population. The ^ 25.3' value is obtained by scaling with 
Aw/viV from the value shown in Figure [TOk for the 10'' Mq 
extended mass. 

We would like to stress that making a definite prediction 
about the A^-dependence of the variance is beyond the scope 
of the current paper. However, we have noted that in both 
cases considered in Figure [TU] the relative variance is of the 
order of unity or larger i.e. the dispersion is of the order of 
the Newtonian peribothron shift. 

The positional uncertainty is currently of the order of 
1 mas. For the highly eccentric orbit of S2 this implies that 
the accuracy with which the peribothron shift can be de- 
tected is of the order of 24'. As can be seen for the case of 
Mi, = 10^ Mq , the shifts are at the limit of the current 
instrumental capabilities if the total enclosed mass was en- 
tirely composed of massive perturbers. The shifts given in 
Equations fF2\ and (I^Bl) can be measured if the accuracy 
is improved by at least one order of magnitude using larger 
telescopes or interferometric methods in the NIR. However, 
considering the variances in the calculated shifts one would 
need to observe more than one stellar orbit in order to infer 
information on the population giving rise to the Newtonian 
peribothron shift. By comparison, the current uncertainty 
in S2's eccentricity is ~ 0.003, and uncertainties in the 
Delaunay angles i and J7 describing its orbital plane are 
~ 50' ([Gillessen et al. 2009b). In both cases, an improve- 
ment of a factor ~ 50 would be required in order to detect 
the changes given in e and d. 

Dynamically-relaxed models of the GC have been criti- 
cized on the grounds that they predict a steeply-rising den- 
sity of old stars insi de ~ 1 pc, while the observations show 
a parsec-scale cor e (jBuchholz et al.l l2009t iDo et al.l l2009t 
iBartko et al.ll2010l ). Dynamically unrelaxed models imply a 



much lower density near Sgr A* and an uncertain fraction 
of stellar-mass black holes (jMerrittI l2010t lAntonini et al.l 
2012). The number of perturbers is so small in these mod- 
els that their effect on the orbital elements of S2 would 
be undetectable for the foreseeable future, barring a lucky 
close encounter with S2. 

In addition to the small amplitude of the perturbations, 
the potential difficulty in constraining N and m comes from 
the nonzero variance of the predicted changes (Figure [TU]). 
The variance in Ao; scales as ^ Aw/-\/iV and would be small 
in the dynamically- relaxed models with N k. lO'^. Another 
source of uncertainty comes from the dependence of the 
amplitude of Aw on 7 (Equation IH]), which is unknown. We 
do not have a good model for predicting the variances in 
|Ae| and AS, but Figure [TU] suggests that the fractional 
variance in these quantities is not a strong function of TV 
or m, and that it is large enough to essentially obscure 
changes due to a factor ~ 5 change in m at fixed M^,. On 
the other hand, considerably more information might be 
available than just Ae and A0 for one star; for instance, 
the full time-dependence of (r, v) for a number of stars. We 
leave a detailed investigation of how well such information 
could constrain the perturber m and iV to a future work. 

4.1. Fighting the limits on the power of stellar orbits 

The results from the previous sub-sections clearly show that 
deriving the net-displacement for an ideal elliptical orbit for 
a single star will not be sufficient to put firm limits on both 
the total amount of extended mass and on the nature of 
the corresponding population. However, the situation may 
be improved if one studies the statistics of the time and 
position dependent deviations along a single star's orbit or 
instead uses the orbits of several stars. 

4.1.1. Improving the single orbit case 

The actual uncertainty in projected right ascension or dec- 
lination, CTp^gj^i^jj, can be thought of as a combination of 

several contributions. Here cr^pparont is the apparent posi- 
tional variation due to the photo-center variations of the 
star while it is moving across the sea of fore- and back- 
ground sources. The scattering process results in a variation 
of positions described by f scattering- Finally, systematic un- 
certainties due to establishing ancf applying an astrometric 
reference frame give a contribution of 0-3^^^^^^^^;^. 

The value of cTpQsJtJon can be measured in comparison 

to the orbital fit. The value of Capparont can be obtained 
experimentally by placing an artificial star into the imag- 
ing frames at positions along the idealized orbit. A reliable 
estimate of cTappaicnt i^ achieved by comparing the known 
positions at which the star has been placed and the posi- 
tions measured in the image frames. As for the case of the 
systematic variations, they can be estimated by investigat- 
ing sources that are significantly brighter or slower than 
the S-stars. Finally, the value that describes the scattering 
process, and therefore gives information on the masses of 
the scattering sources, can be obtained via 



scattering 



position 



apparent 



systematic" 



(31) 



Alternatively, Cg^attering could be measured directly by 
near-infrared interferometry with long baselines. Measuring 



10 



Sabha, Eckart, Merritt, Zamaninasab et al.: The S-star Cluster 



the position of S2 interferometrically as a function of time 
with respect to bright reference objects could allow observ- 
ing the effects of single scattering events. Here the assump- 
tion is that they happen infrequently enough such that one 
can build up sufficient signal to noise on the Cscattcring mea- 
surement provided that the uncertainties in the interfero- 
metric measuring process are sufficiently well known. 

4.1.2. Improving by using several stars 

If scattering events contribute significantly to the uncer- 
tainties in the determination of the orbits, a number of stars 
may help to derive the physical properties of the medium 
through which the stars are moving. While the influence 
of the extended mass imposes a systematic variation of the 
orbits through the Newtonian peribothron shift, the varia- 
tions due to scattering events will be random. This implies 
that for individual stars the effects may partially compen- 
sate or amplify each other. Averaging the results of iV stars, 
that will then essentially sample the shape of the distribu- 
tions shown in Figure El may therefore result in an im- 
provement proportional to N^^/"^ in the determination of 
the extended mass. 

5. Simulating the distribution of fainter stars 

In NSIO we detected three stars that were either previ- 
ously not identified at all (NSl & NS2 stars. Figure 1 in 
NSIO) or only allowed an unsatisfactory identification with 
previo usly known members of t he cluster (S62, as pointed 
out in iDodds-Eden et al.l 120111 ). In addition we have the 
case of the star S3 which was identified in the i^ s-band 
in the e arly epochs 1992 (jEckart fc Genzell I1996D . 1995 
( Ghez et al 1998) and lo st after about 3 ye ars in 1996/7 
( Ghezet al.l[l998h . 1998 (iGenzel et al.ll2000h . We investi- 



gate this phenomenon in our modeling by extrapolating the 
KLF in the inner 1-2 arcsec region, surrounding Sgr A*, to 
stars fainter than the faintest source {K^ — 17.31) we de- 
tected in our 30 August and 23 September 2004 dataset, in 
which Sgr A* shows very low activity (NSIO). In this sec- 
tion we describe the method we use to simulate the distri- 
bution of these faint stars, and the possible false detections 
that can be caused by the combined light of many stars 
appearing in projection to be very close to each other, such 
that they cannot be individually resolved with 8-10 m class 
telescopes. 

The calculations were done by taking all the extra (ex- 
trapolated) faint stars in the Xs-magnitude interval of 18 to 
25. The stars were then distributed in a 23 x 23 grid that 
corresponds to 529 cells. Each cell has the dimensions of 
0.06" X 0.06", i.e. about one angular resolution element in 
i^s-band, this grid, therefore, simulates observations of the 
inner 1.38" x 1.38" projected region surrounding Sgr A*. We 
distributed the faint stars in the grid such that their radial 
profile centered on Sgr A* reproduces that of the stellar 
number density counts of the inner region of the central 
stell ar cluster with a po wer-law index of F = 0.30 ± 0.05 
from lSchodel et all ([20071). This way each cell has a specific 
number of stars that can be inserted into it, with the max- 
imum number of stars being located in the central cell, i.e. 
the peak of the radial profile. Our algorithm fills each cell 
with its specified number of stars by choosing them ran- 
domly from a pool of stars created from the extrapolated 
KLF. The pool is created such that for each i^g-niagnitude 



bin above ~ 18, a number of stars N get their Xg-band 
magnitudes according to the KLF. From this pool of stars 
we then randomly pick objects to fill the cells of the grid 
such that they obey the power-law radial number density 
profile. Then, the fluxes of the stars in each individual cell 
are added up and compared to the value of 0.76 mJy which 
is the flux density of the faintest stellar source in our S-star 
cluster data, i.e. K^ = 17.31 (NSIO). We ran the simulation 
lO'' times in order to get reliable statistical estimates for the 
brightnesses in each resolution cell. Hence we can estimate 
how likely it is to find strong apparent clusterings along 
the line of sight that are brighter than the faintest star we 
identified in the S-cluster (flux larger than 0.76 mJy). 



Table 1. Probabilities of detecting a false star (brighter 
than K^ ^ 17.31) in a 1.38" x 1.38" region. 



Ks-h&uA 


Power-law ir 


dex 


magnitude 








cutoff 


0.19 


0.30 


0.35 




KLF slope = 


0.11 




20.99 


0.0000 


0.1471 


0.1500 


24.67 


0.0285 


0.0292 


0.0224 




KLF slope = 


0.18 




20.99 


0.0848 


0.1286 


0.3016 


24.67 


0.2058 


0.2426 


0.2927 




KLF slope = 


0.25 




20.99 


0.7776 


0.7442 


0.9085 


24.67 


0.9462 


0.9725 


0.9802 



Taking into account the uncertainties of the quantities 
that describe the central S-star cluster we have repeated the 
simulation for a combination of three KLF slopes (0.11, 0.18 
and 0.25), three radial profile power-law indices (F = 0.19, 
0.30 and 0.35) and two i^g-magnitude cutoffs for the extrap- 
olation, 21 and 25 (corresponding to 0.0258 and 0.0009 mJy, 
respectively). Here the brighter cutoff is very close to the 
brightness of the faintest stars that have been detected. The 
choice for the KLF slope satisfies the range of the power-law 
fit F = 0.1 8 lb 0.07. The power- law indices were taken from 
Table 5 of ISchodel et al.l ()2007[ ) for the cusp radial profiles. 

The results of the simulations are summarized in 
Table [T] Three different realizations of a cluster simula- 
tion as well as the average of 10'* simulations are shown 
in Figure [TT] We find that for the measured KLF slope of 
0.18, a measured power-law index of F = 0.3 and a faint Ks- 
magnitude cutoff we obtain a false star in about a quarter 
of all simulations. For steeper KLF and power-law slopes F 
we get this result in more than 70% of all cases independent 
of the cutoff magnitude. 

In Table [D we show the same statistics as in Table [T] 
but for the central cell in the grid, at the projected posi- 
tion of Sgr A*. Also given, in parentheses, is the number 
of stars in the central cell that gives rise to the detection 
of a false star at a distance of less than one angular reso- 
lution element away from the line of sight to Sgr A*. We 
find that for a KLF slope of 0.25 we get a false star in 
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Fig. 11. Upper panel: Three different snapshots of the simulation for the 0.25 KLF slope, power-law index F = 0.3 and 24.67 
it^s-magnitude cutoff. Lower panel: The same as upper but with only the detectable blend stars visible. Right: Average of all the 
10* simulation snapshots for the same setup. 



Table 2. Probabilities of detecting a false star (brighter 
than Ks = 17.31) at the position of Sgr A*. 



ii's-band 

magnitude 

cutoff 



0.19 



Power-law index 



0.30 



0.35 



20.99 
24.67 



KLF slope = 0.11 
0.0000(2) 0.0752(4) 0.0757(4) 
0.0094(6) 0.0103(6) 0.0099(6) 



20.99 
24.67 



KLF slope = 0.18 
0.0423(4) 0.0438(4) 0.1181(5) 
0.0345(15) 0.0591(18) 0.0821(20) 



20.99 
24.67 



KLF slope = 0.25 
0.3120(8) 0.3149(8) 0.5448(10) 
0.3223(59) 0.4756(70) 0.5291(74) 



Notes. The number of stars contributing to the detected flux 
of the false star is given in parentheses for each considered case. 



30% to 50% of aU simulations, independent of the power- 
law index T and the cutoff magnitude. This is consistent 
with the offsets fou nd in different observational ep ochs of 
Sgr A * light curves (jWitz cl et al. 2012; Dodds-Edc n et al.l 
1201 it) . In this case the blend consists of 8 to 74 stars below 
the unresolved background in the S-star cluster region. For 
flatter KLF slopes (i.e. 0.11 and 0.18) we find that a blend 
star only occurs in less than about 10% of all cases, which 
appears to be well below the upper limit found from obser- 
vations. For a KLF slope of a = 0.25 and a number density 
power-law index of F around 0.3 the total number of stars 



in the simulated S-star cluster is a few 1000. This is consis- 
tent with the numb er of main-sequence stars assumed by 
iFreitaeet al.l (I2006D. 



6. Summary and conclusion 

By determining the KLF of the S-star cluster members from 
infrared imaging, using the distribution of the diffuse back- 
ground light and the stellar number density counts, we have 
been able to shed some light on the amount and nature of 
the stellar and dark mass associated with the cluster of high 
velocity S-stars in the immediate vicinity of Sgr A*. 

The amount of light from the fainter S-cluster members 
is below the amount of residual light after removing the 
bright cluster members. One implication could be that both 
the diffuse light and dark mass are overestimated. However, 
while NSIO estimate that only a maximum of one third of 
the diffuse light could be due to residuals from the PSF 
subtraction, we find that faint stars at or beyond the com- 
pleteness limit reached in the KLF can account only for 
about 15% of the background light. Additional light may 
also originate from accretion processes onto a large number 
of 10 Mq black holes that may reside in the central re- 
gion, covered by the S-stars. We find that the stellar mass 
derived from the KLF extrapolation is much smaller than 
the amount of mass that may be present considering the 
uncertainties in the orbital motion of the star S2. Higher 
angular resolution and sensitivity are needed to resolve the 
background light and analyze its origin. 

By investigating the effects of orbital torques due to res- 
onant relaxation, we find that if a significant population of 
10 Mq black holes is present, with enclosed masse s between 
10^ Mq and 10^ Mq (see e.g. lFreitae et al.ll2006[l . then for 
trajectories of S2-like stars, contributions from scattering 
will be important compared to the relativistic or Newtonian 
peribothron shifts. This clearly shows that observing a sin- 
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gle stellar orbit will not be sufficient to put firm limits on 
the total amount of extended mass and on the importance 
of relativistic peribothron shift. In this case only the ob- 
servation of a larger number of stars will allow to sample 
the statistics of the effect, i.e. the distributions in Figure IHl 
However, if the distribution of 10 Mq black holes is cuspy 
then this may become even more difficult (and close en- 
counters should be frequent in this region). 

In general, the inclusion of star-star perturbations al- 
lows us to probe the distribution and composition of mass 
very close to the SMBH simultaneously, if the astromet- 
ric accuracy can be improved by an order of magnitude by 
using either larger telescopes or interferometers in the NIR. 

With measurements and extrapolations of the S-star 
cluster KLF slope, and number density counts with as- 
sumptions on the KLF cutoff magnitude, we can show that 
the contamination for the members of the cluster, and es- 
pecially at the position of Sgr A*, by blend stars is fully 
consistent with measurements. We show that for 8-10 m 
class telescopes the presence and proper motion of faint 
stars close to the confusion limit in the region of the S- 
star cluster is highly contaminated by blend stars. Due to 
the 2-dimensional velocity dispersion of the stars within the 
S-star cluster of about 600 km/s the blend stars will last 
for about 3-4 years before they fade and dissolve. Close to 
the center, we find the probability of detecting blend stars 
at any time is about 30-50%. At the central position the 
change from the appearance of a blend star to the appear- 
ance of another may also give the illusion of high proper 
motions for 8-10 m class telescopes. Such a prime exam- 
ple would be S3, detected close to the position of Sgr A*, 
whic h had both a limited hfetime and high proper mo- 
tion (JEckart fc GenzellfTool iGhez et al.lll998t iGenzel et ahl 
l2000t l. Blending of sources along the line of sight may also 
severely contaminate the proper motion measurements of 
individual stars close to the confusion limit. Only with the 
help of proper motion measurements over time significantly 
longer than 3 years one will be able to derive reliable or- 
bital parameters for a single star. Also, spectroscopy may 
help to resolve blend stars, however, the objects are faint 
and spectroscopy will be difficult. 

These findings clearly demonstrate the necessity of 
higher angular resolution, astrometric accuracy and point 
source sensitivity for future investigations of the S-star clus- 
ter. They would also greatly improve the derivation of the 
amount and the compactness of the central mass as well 
as the determination of relativistic effects in the vicinity of 
Sagittarius A*. 
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